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ABSTRACT 



In this paper we study the production of Li, Be and B nuclei by Galactic cos- 
mic ray spallation processes. We include three kinds of processes: (?) spallation 
by light cosmic rays impinging on interstellar CNO nuclei (direct processes); (ii) 
spallation by CNO cosmic ray nuclei impinging on interstellar p and 4 He (inverse 
processes); and (Hi) a-a fusion reactions. The latter dominate the production 
of 6 ' 7 Li. We calculate production rates for a closed-box Galactic model, verifying 
the quadratic dependence of the Be and B abundances for low values of Z. These 
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are quite general results and are known to disagree with observations. We then 
show that the multi-zone multi-population model we used previously for other as- 
pects of Galactic evolution produces quite good agreement with the linear trend 
observed at low metallicities without fine tuning. We argue that reported dis- 
crepancies between theory and observations do not represent a nucleosynthetic 
problem, but instead are the consequences of inaccurate treatments of Galactic 
evolution. 

Subject headings: Galaxy: evolution — cosmic rays — Galaxy: abundances 

1. Introduction 

The formation of the lightest elements beyond helium ( 6 ' 7 Li, 9 Be, 10,11 B, hereafter LiBeB) 
has persisted as a major problem for nucleosynthesis for many decades. The history is well 
known and we will not review it here (see Spite & Pallavicini 1995; Prantzos, Tosi, & von 
Steiger 1998; Ramaty et al. 1999; Fields & Olive 1999a,b; Vangioni-Flam et al. 2000; Fields 
et al. 2000). Recent observational advances have again brought these elements to center 
stage as data become available for ancient, extremely low metallicity stars in the Galactic 
halo and contemporary, interstellar diffuse gas. 

Several separate issues must be addressed regarding the time development of light el- 
ement abundances in the Galaxy. The Spite plateau appears to be the floor for 7 Li and 
may represent its primordial value. Both isotopes of Li are easily destroyed in stellar inte- 
riors so the likely source for the greater than primordial abundances must be extra-stellar. 
For this reason spallation reactions by Galactic cosmic rays (GCR) interacting with nuclei 
in the diffuse interstellar medium (ISM) have been implicated as avenues for light element 
Galactic synthesis. Recent observations tax the simple spallation models that have p and 
a cosmic rays reacting with stationary interstellar heavy nuclei. As we will show, because 
the heavy nuclei and massive cosmic rays come from the same ultimate source - stars and 
their supernova remnants - the abundance of 6 Li, Be, and B should all scale as Z 2 . This 
prediction is severely at odds with the observational data at low metallicity (Duncan et al. 
1998; Boesgaard et al. 1999). This picture can be modified in several ways. Inverting the 
reasoning, Prantzos et al. (1993) proposed that a change in the slope and lower energy cutoff 
of GCRs over Galactic history can enhance the production of the light nuclei and produce 
an almost linear dependence on Z. Duncan et al. (1997) argued that including fusion and 
inverse reactions, those for which the high energy particles are accelerated CNO spalling 
stationary hydrogen and helium in the ISM, will produce the right metallicity scaling for 
even 7 Li. Fields & Olive (1999a) and Fields et al. (2000) have shown that in a closed box 
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model the early evolution of LiBeB can be explained using the standard rates provided the 
O/Fe ratio is allowed to float with time, increasing at low Z. These calculations were all 
performed assuming only single galactic zones, albeit including simple infall prescriptions. 

Our aim in this paper is to contrast the predictions of two broad classes of schemes for 
treating Galactic chemical and population evolution: simple closed box models, and more 
developed multi-zone, multi-population models. We will show that many of the theoretical 
puzzles noted in the literature are simply paradoxes arising from assumptions in the closed 
one-zone models. These cannot be repaired by more accurate prescriptions of the yields and 
delays due to stellar evolution. They are endemic to the basic feature of the models since the 
system maintains constant mass without feedback and the star formation is a monotonically 
decreasing function of time according to a Schmidt-type power law in gas density (or surface 
density) . 

Once the nuclear reactions involved in the production of LiBeB nuclei are established, 
a Galactic evolution model is needed to compute the ISM abundances evolution and the SN 
rates. Two extreme choices are possible: the simple closed-box model (Tinsley 1980), or a 
more complex nonlinear scheme, that gives satisfactory answers to other aspects of Galactic 
evolution (Ferrini 1991; Ferrini et al. 1992; Pardi et al. 1995). A comparison between the 
two is useful to fix the limits of influence of the hypotheses for the GCR: the observational 
situation of LiBeB reflects simply a nuclear problem of GCR interactions with ISM, and 
hence depends exclusively on knowledge about their interactions. 

2. Cosmic-ray nucleosynthesis 

In this section we consider the nucleosynthesis scenario of spallation reactions of GCR 
on ISM nuclei for the production of LiBeB nuclei. Among all nuclear reactions, there are 
two main production channels of light elements: 

{%) direct processes, i.e. interactions of a p or an a particle of GCR with C,N,0 nuclei 
of the ISM, hereafter represented by p, a + C,N,0 — > LiBeB; 

(ii) inverse processes, i.e. interactions of C,N,0 nuclei of GCR with a p or an a particle 
of the ISM, hereafter represented by C,N,0 + p, a — > LiBeB. 

For Li, fusion reactions must also be considered, i.e. interactions of a particles of GCR 
with a particles of ISM (hereafter represented by a + a — > 6,7 Li). 

The production rate of the light element £ (with £ = LiBeB) at time t, resulting from 
the spallation process i + j — > £, is the convolution of the flux t) of the GCR projectile 
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nuclei i with the numerical abundance of ISM target xf weighted by the spallation cross 



section a^^Ei), 



r id ^ e (t) = xf M (t) / ^(Ei^aij^kiE^Sij^tiE^dEi, (1) 



where is the threshold energy of the reaction. 

We assume that only a fraction Sij->e(Ei) of newly produced light elements is ther- 
malized and incorporated in the ISM, the remaining fraction becoming part of the GCR. 
The reaction products Li, Be and B have, in fact, different fates: they can escape from the 
interaction region if they have a sufficient energy, or suffer a progressive energy loss until 
becoming thermalized. The ability of the ISM to stop the light elements so synthesized de- 
termines the survival probability Sij->e(Ei) of various reaction products against escape from 
the interaction region and for their nuclear destruction that can occur during their slowing 
down. When the incident particles are p or a colliding with heavy nuclei, S^j^^Ej) ~ 1. In 
the opposite case, and for a-a reactions, the trapping fractions are usually calculated in the 
framework of the leaky-box model (see e.g. Cesarsky 1980), and have the form 



where f(Eg[Ei\) is the distribution of energy of the product particle £ (depending on the 
energy Ei of the incident particle), Re(Ee) is the stopping range of the particle £, and A is 
the ISM path length. Here we approximate the trapping fractions as constant, and we adopt 
the results by Mitler (1972), shown in Table 1, that are in agreement with the prescriptions 
by Meneguzzi et al. (1971) and other authors. 

We assume that the GCR flux <&i(E,t) has a separable form, 



where <p(E) is the present-day cosmic ray flux (number of nuclei per unit energy time and 
area, see Sec. 2.1), xf CR (t) is the cosmic ray composition (by number), and is a dimen- 
sionless function that modulates the intensity of the GCR (see Sec. 2.1), with £ = 1 at the 
present epoch. Thus, eq. (1) becomes 





GCR 



(3) 



r id ^e(t) = x { 



GCR 



(t)af M (t)S(t)Vij^eSij^e, 



(4) 



where we define the production coefficients 




(5) 
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The function gives the time dependent intensity of GCR. Following the current 
interpretation, we assume that cosmic rays particles are accelerated by SN remnants, and 
therefore we set 

{«) = y^T. (6) 

V'SN^Gal) 

where ^sn is the global supernova rate (number of supernovae per unit time) and icai — 
13 Gyr is the age of the Galaxy. 

The composition of GCR varies with time according to the evolution of interstellar 
abundances. Thus we set xf CR (t) = x J ^ M (t), where x} SM (t) (i = p, a, C, N, O) is the 
abundance by number of the corresponding element in the ISM, calculated by our model 
of Galactic chemical evolution. Finally, the spallation cross sections (7 i j_ > ^(£'i) for all direct 
and inverse processes are taken from Read & Viola (1984), whereas the cross sections for 
a-a reactions are from Mercer et al. (2001). 



2.1. The energy spectrum of GCR 

A definitive, comprehensive theoretical framework is still lacking that is able to explain 
the origin and mechanism(s) of propagation of Galactic cosmic rays based on the interac- 
tion of relativistic charged particles with the interstellar matter. Approximate semiempirical 
models are widely used to describe the propagation of GCR, often treated as a diffusion pro- 
cess. The "leaky box" model is the most popular: a number of pointlike sources distributed 
throughout the Galaxy emits a flux of fast particles with a broad energy spectrum. The 
particles propagate diffusively within the Galaxy and have a certain probability of escap- 
ing. The homogeneity hypothesis allows a simplified treatment of the corresponding diffusion 
equation for the flux, and the propagation equation can be solved assuming steady state con- 
ditions. The solution depends on several factors: the energy losses by ionization, the mean 
free path for escape from the Galaxy, the mean free path for nuclear reactions, the assumed 
spectrum at the source, etc. Energy loss resulting from continuous ionization accounts also 
for the probability that a particle from the source or produced by nuclear reactions during 
the propagation of GCR can survive to become part of the flux. 

In this work we adopt the propagated demodulated spectrum obtained by Ip & Axford 
(1985), which we approximate as 

<p(E) = A(E + E )-\ (7) 

where E is in MeV, E = 750 MeV, and 7 = 2.7. For E > E this spectrum becomes a 
simple power law in energy that allows to adopt the absolute flux normalization A = 3.88 x 
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10 6 cm- 2 s- 1 MeV- 1 obtained by Wiebel-Sooth & Biermann (1999) from data at E = 1 TeV. 
Alternative choices are available in the literature (Gloecker & Jokipii 1967; Goldstein et 
al. 1970; Dodds et al. 1976; Reeves & Meyer 1978), but their agreement with this assumed 
spectrum is good. Our adopted spectrum also agrees well with the more recent results quoted 
by Webber (1998); discrepancies of around 50% are found below 100 MeV where solar wind 
modulation is important, but above that the general agreement is within 10%. We do not 
include an extended low energy tail that has been invoked in previous models, but as we will 
show in section 4, the full treatment for evolution does not require this. The values of Vi^i 
for the reactions considered, evaluated with the GCR spectrum and nuclear cross section 
just described, are shown in Table 2. These production coefficients contain all the necessary 
nucleosynthesis information. 



3. Analytic Results Models for Light Element Evolution in a Closed Box 

Formalism 

The simplest representation of Galactic evolution is the single zone closed box scheme, 
also called the Simple Model (SM, see e.g. Tinsley 1980). It provides a useful, although 
limited, test bed for model assumptions and here we use it as a guide to more complete 
chemical evolution modeling that forms the next section. The main advantage of this closed 
box picture is that we can derive analytical expressions for the Li, Be and B abundances, 
evaluate orders of magnitude of the most significant quantities, and gain an idea of their 
trends. It should be kept in mind, however, that the model is severely restricted: there can 
be no net flow of mass through the zone, nor is there coupling with any other parts of the 
Galaxy. 

In the SM, the Galaxy is a spatially homogeneous, isolated closed volume of fixed total 
mass M . The gas mass M g decreases with time from star formation at a rate ip, 

where R is the average mass fraction, weighted by the initial mass function (IMF) <p(m), 
returned to the ISM by stars at the end of their evolution. In the following, we employ the 
instantaneous recycling approximation (IRA), and assume that stars process interstellar gas 
and replenish it on a timescale shorter than the typical large scale dynamical timescales. 
The rate of metal production by stars is: 



(9) 
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where Z is the metallicity and Pz is the mass fraction of heavy elements newly produced and 
ejected by a stellar generation weighted by the IMF. As usual in such treatments, we make 
no distinction between the stellar and interstellar abundances. Here we adopt the numerical 
values obtained by Galli et al. (1995) from the same stellar and nucleosynthetic data adopted 
in the evolution model of Ferrini et al. (1992) R = 0.21 and Pz = 7.9 x 10~ 3 . In keeping 
with the restrictions of the IRA, massive stars are the principal nucleosynthetic agents and 
since the IMF remains unaltered with time, both R and Pz are held constant. 

In a closed system, since the total mass M must remain fixed, M g and t are conjugate 
variables and eqs. (8) and (9) yield 



z -rs'"Uj- (10) 

with Z(0) = 0. Notice that we can include the effects of global mass loss through star- 
induced processes, i.e. winds, supernovae, etc., by inserting here a term proportional to ip. 
The net effect is the same as reducing R. On the other hand, gas infall is not so easily 
treated. It actually violates the basic assumption of the closed box: there is a source of 
matter from outside the system. Thus M g and Z are explicit functions of time and must be 
found numerically. 

Light elements are mainly destroyed in stars during their evolution, except for limited 
sources of specific isotopes (during the AGB stage and possibly from some novae). On the 
other side of the coin, they are produced in the ISM at a rate ri(t) = £\ . rjj_^(£), therefore 
the abundance (by number) xg of the element £ is determined by the equation 

j t (x e M g ) = -x^ + r e M g , (11) 

or, using equation (8), by the equation 

M g ^- = -x e R^ + r e M g . (12) 

To solve equation (12), it is necessary to evaluate the term r^(t), using equation (4). 
Following convention we assume that the normalized supernova rate is proportional to the 
SFR, £(£) = V'SM^V'^'SN^Gai) — V'(*)/V'(*Gai)) making no distinction between supernovae 
types I and II. This is a good approximation because the stars that end their life as SN 
II are high mass stars and then characterized by short lifetimes with respect to the age of 
the Galaxy. We only use supernovae to scale the production rate for cosmic rays, not for 
the overall metallicity of the gas (for which the two types are implicitly combined in the 
IRA coefficients). As we will describe in the next section, the distinction between SN types 
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is made explicitly in the fully nonlinear chemical evolution models. With our assumption 



x 



GCR 



X 



ISM 



(t), the production coefficients r> for both direct and inverse processes scale 



linearly with the gas metallicity Z, since the abundance of p and a remains nearly constant 
with time. For fusion ol—ol reactions the production coefficient is practically independent on 
Z. Therefore, for direct p and a induced reactions and inverse CNO induced reactions, we 
obtain 

z (t) _ \ ^(*Gai)*o 



-re(t Ga \) = 



-Mt)Z{t), 



(13) 



Z(t Gal ) ,K " y M Z(t Gal ) 
where we have defined a timescale t = M /ijj(t Ga ]). For 6 Li and 7 Li we must also include 
fusion ol-ol reactions, 

(14) 



M 

Although the same analytic expression describes both direct and inverse processes, keep in 
mind that under our assumptions the Z dependence in the first case represents the instan- 
taneous interstellar abundances, while in the second case it is the GCR composition. 

With these definitions, it is straightforward to solve equation (12) together with equation 
(8) and (9): 



+ 



: x e ,BB exp 
re(t G3l )t P z 



exp 



(1 - 2RfZ{t GaX ) 
for direct and inverse processes, and 



Pz 



1 + 



1-2R 



exp 



1-R 



x e (Z) 



,bb exp 



Pz 



+ 



r t (t 



GaXjto 



1 - 2R 



exp 



Pz 



exp 



1-R 



-Z 



(15) 



(16) 



for fusion reactions. Any contribution from Big Bang nucleosynthesis is included as an 
initial abundance x^bb- Only for 7 Li the abundance depends on all of the processes we have 
discussed, including a significant cosmological contribution. 

Equations (15) and (16) show that the light elements evolution is determined by the 
combination of stellar destruction, which reduces the initial Big Bang abundance, and cosmic- 
ray production that is proportional to r^(t Gal )t - Thus, according to the SM, all the physics of 
Galactic evolution is parameterized by the single quantity to while the remaining quantities 
in equations (15) and (16) depend only on stellar properties. 

In the absence of infall, the star formation rate must exponentially decrease with time 
since inert remnants inevitably result from stellar evolution. With this assumed rate, it is 
straightforward to solve eq. (8) and (9), obtaining 

M g (t) = M e~ t/T , (17) 
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and 

m = (T^ 4 ' (19 > 

From eq. (19) we obtain r ~ 4.2 Gyr by requiring that by t — t Q — 8.5 Gyr the metallicity 
has reached the solar value Z Q = 0.02. For tQai — 13 Gyr, eq. (19) also gives Z(tQ a \) ~ 1.5Z Q . 
Notice that Sandage (1986) estimated r ~ 4 Gyr. Then for to, we obtain 

t = J}p— = (1 - i?)re^ l/r ~ 72 Gyr. (20) 

The coefficients r^(t Ga i) are given by eq. (4). The production coefficients Vij->e and the 
trapping fractions Sij^g are given by Table 1 and 2, respectively. For both GCR and ISM 
at t — ioai we have adopted a solar composition, xh = 0.91, :thc = 0.089, xq = 3.30 x 10 -4 , 
xn = 1-02 x 10~ 4 , and xo = 7.73 x 10~ 4 (Grevesse & Sauval 1998). The values of r^tcai) 
obtained in this way are listed in Table 3. With these, and the value of to given by eq. (20), 
we obtain the curves shown in Figures 1-4. 

For Li evolution, shown in Figures 1 and 2, the predicted abundance remains very close 
to the primordial value, represented by the Spite plateau, up to [Fe/H] ~ —1.4. At higher 
metallicities the predicted abundance increases about tenfold, although remaining signifi- 
cantly below the upper envelope of the observational data and the solar system meteoritic 
value. We will see in Section 4 that a more realistic model of Galactic evolution considerably 
reduces the predicted Li abundance at high metallicity, thus highlighting the need to include 
an additional (stellar) source of Li at later Galactic epochs. In Fig. 1, we see that the main 
contribution to the rise at [Fe/H]> —2 comes from fusion reactions that depend linearly 
on the relative fraction of GCR a particles. If, as reported in e.g. Meyer et al. (1998) and 
Wiebel-Sooth & Biermann (1999), these are depleted relative to CNO by as much as a factor 
of 5 without a change in spectrum, then the rise is delayed and the enrichment of Li reduced 
as displayed. The SM predictions for Be and B (Figures 3 and 4) produce reasonable fits to 
the observed abundances near Z ~ Z Q , but fail at low metallicity. This is simply seen from 
equation (15) which, for Z < Z , gives the approximate behavior of the abundance of Be 
and B at low metallicity: 

MZ) - if^^V, (2D 

showing a quadratic dependence on Z, at variance with the nearly linear trend up to [Fe/H] 
~ —1.3 (see Boesgaard et al. 1999) in the observational data. In contrast, Li (= 6 Li+ 7 Li) is 
predicted to have a linear dependence with Z, 

xu ^ a; L i,BB H 5 Z (22) 
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or, 




(23) 



where 



k ~ 



xkPz 



(24) 



and r Li (t Gal ) = r 6 (t Gal ) + r 7 (t Gal ) ~ 9.8 x 10~ 28 s _1 for fusion reactions (see Table 3). In 
eq. (24) we have neglected the small effect of astration of primordial Li represented by the 
term — Rx^^b in eq. (23). The value of k has been determined observationally by Ryan et al. 
(2000) in a relatively small sample of carefully selected low metallicity stars (Ryan, Norris 
& Beers 1999). For a parametric linear fit they find k — (1.1 ± 0.7) x 10~ 8 . With x H = 0.91 
the SM with our parameters predicts k ~ 6.2 x 10~ 9 . This accordance indicates that the 
spallation scenario is likely correct, not that the closed box model is necessarily appropriate, 
because this expansion is correct only for near solar abundances and, as we have seen, the SM 
predictions are not consistent with each other. We therefore pass to the more complex model 
and show how the light element evolution can be better modeled, though at the expense of 
a considerably more complex treatment for the Galaxy. 

4. Results for Light Element Evolution in a Multizone-Multi Population 



We now consider the production of LiBeB nuclei in the framework of a more complex 
model of Galactic evolution. We recall briefly some important features of the model adopted 
(for a more detailed discussion see Ferrini et al. 1992, 1994). This multiphase and multi-zone 



model is able to follow the abundances of H, D, 3 He, 4 He, 12 C, 13 C, 14 N, 16 0, 20 Ne, 24 Mg, 28 Si, 



32 S, 40 Ca, and 56 Fe and has recently been extended to s- and r-process elements (Travaglio 
et al. 1999). The model follows the interconnected evolution of three different regions in the 
Galaxy: halo, thick disk, and thin disk (multizone treatment). Relaxing the IRA takes into 
account detailed stellar nucleosynthesis as well as finite stellar lifetimes. The multiphase 
approach considers the different phases of the matter in the Galactic system: (i) stellar 
populations; (ii) stellar remnants; (Hi) interstellar matter phases. A fundamental feature of 
the model is self-regulation through the internal processes that allow for the transformation 
of matter from one phase to another. These interactions produce the time dependence of 
the total mass fraction in each phase and the chemical abundances in the ISM and in stars. 
The star formation rate ip(t) is determined self-consistently by the interaction among the 
different phases of matter. This is very important because in the majority of the existing 
models, such as the SM described in Sec. 3, the formal time dependence for the SFR is 
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assumed a priori The initial mass function adopted in this model is based on the analysis 
of Ferrini (1991) of the fragmentation of molecular clouds in the solar neighborhood, and is 
constant in space and time. 

We have included the nucleosynthesis processes responsible for the evolution of LiBeB 
nuclei in this general Galactic model but we emphasize that no other modification has been 
made to the model parameters or equations. The results we report here require the input 
from our standard version that has been successfully used to analyze a number of Galactic 
evolution problems (see e.g. Shore & Ferrini 1995). The nonlinear evolution equations that 
determine the Li, Be and B abundances xi (by number), including spallation and fusion 
process, are: 

—7— = r ijH - x e ,H : , 25 

at g H + c H 

dxer Wt , \ 9h 

—rr- = r i:T - x e , T : + jH\xt,H - xi,t) : , (26) 

at 9t + ct 9t + ct 

dx e,D W D g T 

— — = r e , D - x e ,D : + jT{Xe,T - x e>D ) ■ , (27) 

at g D + c D g D + c D 

where £ = 6 Li, 7 Li, 9 Be, 1X B is the light element under consideration. Here /# and fx are, 
respectively, the infall rates for gas flowing from the halo into the thick disk, and from 
thick disk into the thin disk: /# = 0.10 and fx = 0.0065 are the choices made for the 
solar neighborhood, that reproduce the mass ratio of the visible mass in the three zones. 
The quantities gi, and cl are the gas and cloud fractions in the halo, thick disk, thin disk, 
respectively, where L = H, T, or D denotes the Galactic zone. The destruction terms in the 
three zones are represented by W L (t) = ^Wj^^), where 

W itL (t) = J2 QvMXjlt - T Z (m)\i\) L [t - T Z (m)}dm, (28) 

is the mass fraction injected in the ISM as element i from stars ending their evolution at the 
time t. In order to compute the nucleosynthesis production of a star of mass m, the astration 
matrix Qij(m) has been introduced, following Talbot & Arnett (1973). The {ij} element of 
this matrix represents the fraction of the star mass initially in the form of chemical species 
j, transformed and ejected as chemical species i, Qij(m) = m^ p /rrij. The quantity Qij(m) is 
obtained weighting the contribution of the normal evolving stars, including stars exploding 
as SN II and the binary system ending as SN I. As we discussed in the previous section, 
for the closed box the two types have been combined for their contributions to cosmic ray 
production but the full chemical evolution model accounts for the distinctions between these 
two types of supernovae. 
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In Figures 2-4 we show the abundances predicted by our model, obtained following the 
nucleosynthesis prescriptions described in Section 2. The different trends for the metallicity 
as function of time contribute as well to the differences apparent as shown in Figure 5. In 
the multiphase model, Z is not a linear function of time. In fact, the IRA is relaxed and 
the consequences are particularly notable at low Z, when the finite stellar lifetime is most 
important. After t ~ 2 Gyr, (or logZ ~ —2.5) the SM provides a good approximation of 
the numerical model. This is just the range, [Fe/H]> —0.8, where the abundances of LiBeB 
elements predicted by the two schemes are most similar. The displacement is essentially in 
[Fe/H] for this reason. 

Another important question for chemical evolution models is the normalization in the 
predicted abundances of LiBeB elements. It is common practice (see for example Prantzos 
et al. 1993; Fields & Olive 1999a,b; Vangioni-Flam & Casse 1999) to normalize the final 
abundances of GCR-only isotopes 6 Li, 9 Be, 10 B to the observed abundances at [Fe/H] = 0. 
Consequently the GCR yields of 7 Li and n B are scaled by a factor that is the average of the 
three scaling factors of these isotopes. We stress that the results of our models, shown in 
Figures 2-4, are obtained without any ad hoc normalization. The most evident and relevant 
feature of the new curves is its improved agreement with the observations: whereas the 
simple model yields a quadratic dependence of [Be/H] and [B/H] with [Fe/H], the multiphase 
Galactic evolution model more nearly approaches the observed approximately linear trend 
for increasing values of metallicity. Thus, the slope for the abundance of Be and B vs. 
[Fe/H] is ~ 2 for the halo, ~ 1.4 for the thick disk, ~ 0.7 for the thin disk. For comparison, 
Boesgaard et al. (1999) determined empirically for Be vs. [Fe/H] a best-fit slope ~ 1-1.3 for 
-3.0 < [Fe/H] < -1, and ~ 0.6-0.7 for -1.0 < [Fe/H] < 0.1. It worth stressing this result 
because it is different from the conclusions obtained previously by other authors (see for 
example Vangioni-Flam & Casse 1999; Fields & Olive 1999a). Even treating only the direct 
processes (i.e. the "GCR standard" as named in Vangioni-Flam & Casse 1999), it is possible 
to reproduce quite well the observed trend of Be and B with [Fe/H]. From this trial we can 
conclude that in our scenario of Galactic evolution the "GCR standard" does not lead to 
a quadratic dependence of Be and B vs Fe, in contrast to, for example, Vangioni-Flam & 
Casse (1999) and Fields & Olive (1999a). 

The model reaches the solar photospheric B abundance log(B/H)= —9.45 (Grevesse & 
Sauval 1998) at [Fe/H] = 0, whereas it exceeds the solar Be abundance log(Be/H)= —10.85 
(Chmielewski, Miiller, & Brault 1975). This trend also appears in other published results 
(see e.g. Lemoine et al. 1998) that overproduce Be relative to its meteoritic and photospheric 
values. This discrepancy with the observed Be and B abundances at higher metallicity may 
be partially due to stellar depletion. The systematic departures at more recent times of 
the numerical predictions from the observations suggests this interpretation. In Table 4 we 
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compare the model results at solar metallicity with the solar system data (Hanon et al. 1999). 
We remind the reader that the model results refer to ISM abundances, whereas spectroscopic 
observations can only provide stellar photospheric values. A comparison of the theoretical 
curves and the Be abundance at the time of formation to the Solar System (meteoritic value) 
and in the ISM would represent a more appropriate test of the model. Although the ISM 
value 4.6 Gyr ago is unknown, Howk et al. (2000) provide a lower limit on the contemporary 
interstellar logB/H > —9.60. If we assume that this is the same difference as between the 
models and the stellar observations, there are two possible explanations. The more radical 
is that there could be indications of a mixing mechanism - presently unexplained - for stars 
younger than about 1 to 2 Gyr. This is based on the metallicity at which the first systematic 
deviations appear in the Be and B evolution models. We re-iterate that he models have not 
been normalized nor have any parameters used for the stellar population and heavy metal 
evolution been altered to obtain agreements between the observed and predicted light element 
abundances. We favor a more prosaic explanation: that the Be and B abundances have been 
systematically underestimated for the higher metallicity stars, likely due the combined effects 
of NLTE and the "missing opacity" problem most recently outlined by Bell et al. (2001) and 
references therein. If so, assuming the appropriateness of our embedded assumptions about 
the production rates, the models point to the systematic correction that is needed and could 
also serve as a constraint in the search for the responsible opacity sources without the need 
to invoke new production sources for the light elements at late times in galactic history. 

Direct processes and the fusion reaction do not reproduce the observed values for 7 Li 
at high [Fe/H]. This is similar to the results of the SM and also to those obtained in the 
literature. It is well known that the case of 7 Li is more complex: lithium is not a GCR- 
only element and stellar production cannot be ignored (see e.g. Travaglio et al. 2001). For 
11 B/ 10 B, the models predict an essentially constant ratio, 2.4, from spallation processes alone. 
A possible solution is the production of 11 B by //-induced spallation of 12 C in SN II (Woosley 
& Weaver 1995) but this and other likely contributing sources have not been included in the 
calculation. From Figures 2-4, we note a decrease of the slope of the the halo and thick 
disk curves at low [Fe/H] in the cases of Be and B. This is in better agreement with the 
observational data in this range than the results obtained including only direct processes. 
The contribution of inverse processes is more significant at low metallicity, because during 
the first evolutionary phases of the Galaxy the abundance of interstellar CNO target nuclei 
is very low. For Li, the contribution of inverse processes is marginal with respect to the 
pre-Galactic Big Bang production. 
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5. Discussion and Conclusions 

The SM is a closed box model with IRA that provides a useful test for more complex 
models of Galactic evolution because its predicted metallicity is the highest possible. This is 
due to the approximate treatment (IRA) of gas return, which maximizes metal production 
from stars and the efficiency of its re-incorporation into stars, and to the absence of gas 
dilution by infall. Therefore, in the usual abundance vs. [Fe/H] plots, SM predictions 
roughly determine the right boundary of the region occupied by more complex models for 
which these assumptions are relaxed; our multiphase model of Galactic evolution relaxes 
the IRA approximation and allows infall of gas allowing us to follow simultaneously the 
evolution of three Galactic zones (halo, thick disk and thin disk). This fact is evident from 
the comparison of the results of the SM with those of our model shown in Figures 2, 3 and 
4. For the halo component, the two models give similar results. This is because the halo 
behaves like a closed box, even though it loses gas to form the thick disk. Here the loss of 
gas is not what matters but the fact that there is no feedback on the halo from the other 
Galactic zones. Thus, the LiBeB isotope productions in the halo follow the prediction of 
the SM but are shifted to lowest values of [Fe/H] because the IRA overestimates the actual 
metal production. The slope of log(Be/H) and log(B/H) vs. [Fe/H] for the halo is close to 2, 
the value for the SM. Indeed, this quadratic dependence of the abundance of these elements 
on metallicity in the halo is a natural consequence of our halo picture. 

On the other hand, for the thick and thin disk, the IRA is not the only difference with 
the SM and the couplings halo-thick disk, and thick disk-thin disk play important roles 
which are, of course, ignored in closed box models. Thus, in addition to an overall rightward 
shift of the abundances vs. [Fe/H] with respect to the SM the fact that the thick disk 
receives chemically enriched gas from the halo, and the thick disk from the thick disk, allows 
progressively shallower slopes of the abundance curves vs. [Fe/H]. Our goal here is not much 
to discriminate among different stellar populations in the Galaxy. Rather, we deliberately 
loosely use the terms "halo" , "thick disk" and "thin disk" to describe coupled Galactic zones 
forming in sequential cascade, and we emphasize the different behavior of such a complex 
system with respect to a homogeneous one- zone model, no matter how sophisticated. 

An additional complication in the interpretation of the results is the customary (but 
unavoidable) way of plotting the abundance of chemical species not as function of time, as 
they are calculated, but as function of metallicity Z or [Fe/H], to allow a comparison with 
spectroscopic data. Whereas time and metallicity scale linearly in the case of the SM and the 
use of either variable is equivalent, the same is not true for our multiphase model. First, the 
relation of metallicity with time is not linear. Second, each Galactic zone follows a different 
enrichment history. For example, during the first 3 Gyr of Galactic life (where [Fe/H] rises 
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from minus infinity up to about —0.7) the abundance of LiBeB elements is higher in the 
halo, followed by the thick disk and the thin disk. This simply reflects the behavior of the 
star formation rates in the three zones (see Fig. 5 of Pardi et al. 1995); for the same reason, 
[Fe/H] follows the same ordering (see also their Fig. 6), so that a given value of [Fe/H] is 
reached first by the halo and later on by the thick disk and the thin disk. Since LiBeB 
abundances grow more rapidly than Fe with time, however, their ordering in the abundance 
vs. [Fe/H] plot is reversed. At a specified value of [Fe/H], the halo has the lowest abundance, 
and the thin disk has the highest. These facts should be kept in mind when interpreting our 
results for LiBeB shown in Fig. 2-4. 

In conclusion, we show that the multi-zone multi-population model we used previously 
for other aspects of Galactic evolution produces quite good agreement with the linear trend 
observed at low metallicities as consequences of multi-zone coupling of a galactic chemical 
evolution model and detailed accounting for the feedback between phases without the instan- 
taneous recycling approximation or time dependent variations in the cosmic-ray spectrum 
or the initial mass function. No special mechanisms are needed to produce the increase in 
abundances at late times in Galactic history (e.g. novae producing excess Li, variations in 
the low energy part of the cosmic ray spectrum, etc.). Indeed, the overproduction of the 
light elements in our calculations suggests that the abundances of these elements derived 
from stellar photospheric observations may be underestimates. We argue that reported dis- 
crepancies between theory and observations do not represent a nucleosynthetic problem, but 
instead are the consequences of inaccurate treatments of Galactic evolution. 
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Table 1. Trapping fractions S id ^ e 



i j £ Sij^e 



p C,N,0 6J Li, 9 Be, 10 ' n B 1 

a C,N,0 6 - 7 Li, 9 Be, 10 - n B 0.86 

a a 6 - 7 Li 0.46 

C,N,0 p,a 6 ' 7 Li, 9 Be, 10 > n B 0.25 



Table 2. PRODUCTION COEFFICIENTS Vij^e (IN 10~ 25 S _1 ) 



i,j or j,i 


6 Li 


7 Li 


9 Be 


io B 


n B 


p,C 


2.64 


4.82 


0.723 


5.00 


13.3 


p,N 


4.21 


2.52 


1.08 


2.94 


5.93 


p,0 


3.03 


4.75 


0.957 


3.15 


6.91 


a,C 


11.5 


15.2 


3.84 


13.0 


24.1 


a,N 


4.17 


4.80 


1.83 


9.26 


15.2 


a,0 


3.93 


5.11 


1.86 


0.87 


11.2 


cc,ai 


0.74 


1.94 
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Table 3. Production rates rjj_^(£ G ai) for the closed-box model (in 10 30 s 1 ) 



i,j 


6 Li 


7 Li 


9 Be 


io B 


n B 


direct processes 












p,C 


79.3 


145 


21.7 


150 


399 


p,N 


39.0 


23.3 


10.0 


27.2 


55.0 


p,0 


213 


334 


67.3 


222 


486 


total p,CNO 


331 


502 


99.0 


399 


940 


a,C 


29.0 


38.4 


9.70 


32.8 


60.9 


a,N 


3.26 


3.75 


1.43 


7.23 


11.9 


a,0 


23.2 


30.2 


11.0 


5.15 


66.3 


total a,CNO 


55.5 


72.4 


22.1 


45.2 


139 


a, a 


270 


707 








inverse processes 












C,p 


19.8 


36.2 


5.43 


37.5 


99.9 


N,p 


9.77 


5.85 


2.51 


6.82 


13.8 


0,p 


53.3 


83.5 


16.8 


55.4 


122 


total CNO,p 


82.9 


126 


24.7 


99.7 


236 


C,a 


8.44 


11.2 


2.82 


9.54 


17.7 


N,ct 


0.94 


1.09 


0.42 


2.10 


3.45 


0,a 


6.76 


8.79 


3.20 


1.50 


19.3 


total CNO,a 


16.1 


21.1 


6.44 


13.1 


40.5 
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Table 4. Solar System and ISM LiBeB Values Compared to Model 



Value ( 6 Li+ 7 Li)/H 9 Be/H ( 10 B + n B)/H 



Chondrites ( a ) 


5.1 x 


io- 


-10 


2.1 


X 


10" 


-11 


3.3 x 


10" 


-10 


Orgueil ( a ) 


2.0 x 


10" 


-9 


2.6 


x 


10" 


-11 


7.6 x 


10" 


-10 


Solar Photosphere ( b 


) 1.3 x 


io- 


-11 


1.4 


x 


10" 


-11 


3.5 x 


10" 


-10 


ISM 


> 3.7 


x 10- 10 ( c 


) ... 








> 2.5 


x 10- 10 ( d ) 


Model: [Fe/H]=0 


1.1 x 


10- 


-9 


6.0 


x 


10- 


-11 


5.3 x 


10- 


-10 



References. - ( a ) Hanon et al. (1999), ( b ) Grevesse & Sauval (1998), 
( c ) Lemoine et al. (1993), ( d ) Howk et al. (2000) 



-22 - 




Fig. 1. — Contributing physical processes to Li evolution (= 6 Li+ 7 Li) according to the closed 
box model. Dotted line: Big Bang nucleosynthesis and stellar astration; short- dashed line: 
fusion a-a reactions; long-dashed line: direct processes; dash-dotted line: inverse processes. 
The upper most curve shows the result obtained with a standard (solar) GCR composition 
as described in Sect. 3, the lower curve employed a reduced a contribution in the GCR flux 
by a factor of 5 (see text). 
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Fig. 2. — Evolution of Li (= 6 Li+ 7 Li) according to the closed box model (thin dotted line) 
and our numerical model (thick short-dashed line, halo; thick long-dashed line, thick disk; 
thick solid line, disk) including the contribution of fusion, direct and inverse reactions. The 
meteoritic and the solar photospheric Li abundances (from Grevesse & Sauval 1998) are 
indicated by an open circle and an asterisk, respectively. The lower limit on the ISM abun- 
dance of Li determined by Lemoine et al. (1993) has been placed at [Fe/H]= 0.2. Stellar 
photospheric data are from Bonifacio & Molaro (1997) (open squares), Balachandran (1990) 
(filled squares), Fulbright (2000) (filled circles), and Chen et al. (2001) (filled triangles). 
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[Fe/H] 



Fig. 3. — Evolution of 9 Be according to the closed box model (thin dotted line) and our 
numerical model (thick short-dashed line, halo; thick long-dashed line, thick disk; thick solid 
line, disk) including the contribution of direct and inverse reactions. The meteoritic (from 
Grevesse & Sauval 1998) and the solar photospheric (from Chmielewski, Miiller, & Brault 
1975) 9 Be abundances are indicated by an open circle and an asterisk, respectively. Stellar 
photospheric data are from Boesgaard et al. (1999) (filled triangles), Primas et al. (2000a) 
(G64-12, open square), and Primas et al. (2000b) (LP815-43, filled square). 
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[Fe/H] 



Fig. 4. — Evolution of B (= 10 B+ 11 B) according to the closed box model (thin dotted line) and 
our numerical model (thick short-dashed line, halo; thick long-dashed line, thick disk; thick 
solid line, disk) including the contribution of direct and inverse reactions. The meteoritic 
and the solar photospheric B abundances (from Grevesse & Sauval 1998) are indicated by 
an open circle and an asterisk, respectively. The lower limit on the ISM abundance of B 
determined by Howk et al. (2000) has been placed at [Fe/H]= 0.2. Stellar photospheric data 
are from Cunha et al. (2000) (filled squares), Primas et al. (1999) (open triangles), Duncan 
et al. (1997) (filled circles), and Garcia Lopez et al. (1998) (empty squares). 
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Fig. 5. — Comparison of the metallicity Z predicted by the SM (dashed line) and by our 
model (solid line) as function of time. The solar symbol indicates the solar metallicity at 
the epoch of formation of the Sun, t = 8.5 Gyr, assumed in our model. 



